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There is longstanding fundamental interest in 6-fold coordinated dS (t%g) transition metal com¬ 
plexes such as [Ru(bpy) 3 ]^^ and Ir(ppy) 3 , particularly their phosphorescence. This interest has 
increased with the growing realisation that many of these complexes have potential uses in applica¬ 
tions including photovoltaics, imaging, sensing, and light-emitting diodes. In order to design new 
complexes with properties tailored for specific applications a detailed understanding of the low- 
energy excited states, particularly the lowest energy triplet state, Ti, is required. Here we describe 
a model of pseudo-octahedral complexes based on a pseudo-angular momentum representation and 
show that the predictions of this model are in excellent agreement with experiment - even when the 
deviations from octahedral symmetry are large. This model gives a natural explanation of zero-field 
splitting of Ti and of the relative radiative rates of the three sublevels in terms of the conservation of 
time-reversal parity and total angular momentum modulo two. We show that the broad parameter 
regime consistent with the experimental data implies significant localization of the excited state. 


I. INTRODUCTION 

Six-fold coordinated (t^g) transition metal com¬ 
plexes, such as those shown in Fig. ,b, share many 
common properties. These include their marked similar¬ 
ities in their low-energy spectra [T] , cf. Table and the 
competition between localization and delocalizsation in 
their excited states [2]. Beyond their intrinsic scientific 
interest, understanding and controlling this phenomenol¬ 
ogy is further motivated by the potential for the use 
of such complexes in diverse applications including dye- 
sensitized solar cells, non-linear optics, photocatalysis, 
biological imaging, chemical and biological sensing, pho¬ 
todynamic therapy, light-emitting electro-chemical cells 
and organic light emitting diodes [nisHz]. As many of 
these applications make use of the excited state prop¬ 
erties of these complexes a deep understanding of the 
low-energy excited states, particularly the lowest energy 
triplet state, Ti, is required to enable the rational design 
of new complexes. 

Coordination complexes where there is strong spin- 
orbit coupling (SOC) present a particular challenge to 
theory because of the need to describe both the ligand 
field and the relativistic effects correctly. There has 
been significant progress in applying relativistic time- 
dependent density functional theory (TDDFT) to such 
complexes; but significant challenges remain, for example 
correctly describing the zero-field splitting [SUIT]. There 
has been less recent focus on the use of semi-empirical ap¬ 
proaches, such as ligand field theory [HlIIlHIl]. However, 
semi-empirical approaches have an important role to play 
m- Firstly, they provide a general framework to un¬ 
derstand experimental and computational results across 
whole classes of complexes. Secondly, when properly pa- 
rameterised they can provide accuracy that is compet¬ 
itive with first principles methods. Thirdly, they can 
provide general design rules that allow one to effectively 
target new complexes for specific applications. 


A long standing question in these complexes is whether 
the excited state is localized to a single ligand or delo¬ 
calized [2] . The main semi-empricial approach to under¬ 
standing organometallic complexes is ligand field theory. 
Once all of the spatial symmetries are broken there is lig¬ 
and field theory is limited to a perturbative regime near 
approximate symmetries, this makes an accurate descrip¬ 
tion of localised excited states challenging. 

In this paper we describe a semi-empirical approach, 
based on the pseudo-angular momentum approach that 
has found widespread use in, e.g., interpreting electron 
paramagnetic resonance experiments. We derive con¬ 
servation laws based on the total angular momentum 
(pseudo plus spin) that apply even when the pseudo- 
octahedral and trigonal symmetries are strongly broken. 
These conservation laws imply selection rules for radia¬ 
tive emission. We show that this model reproduces the 
experimentally measured trends in the radiative decay 
rates and excitation energies for all of the complexes for 
which we have data to compare with in the literature. 
These trends are insensitive to the parameters of the 
model studied. Finally, we show that for the wide pa¬ 
rameter range compatible with experiment the pseudo- 
angular momentum model predicts significant localiza¬ 
tion of the excited state. 


II. THE PSEUDO-ANGULAR MOMENTUM 
MODEL. 

It has long been understood [16] that the three-fold 
degenerate states can be represented by an / = 1 pseudo- 
angular momentum. Perhaps the best known example of 
this are the t 2 g states of a transition metal in an octa¬ 
hedral ligand field. In the complexes considered here 
the t 2 g orbitals are filled, whereas the e^-orbitals are high 
lying virtual states. Therefore we only include the t 2 g or¬ 
bitals in the model described below. 

The complexes listed in Table [I| have 6-fold coordi- 
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FIG. 1: The structures of two important pseudo-octahedral 
transition metal complexes: a) [Ru(bpy) 3 ]^^ and b) Ir(ppy) 3 , 
where bpy is bipyridine and ppy is 2-phenylpyridyL Sketches 
of the c) TT and d) tt* orbitals of a bpy ligand with the reflec¬ 
tion plane marked by the dashed line. It is clear that these 
correspond to the bonding and antibonding combinations of 
singly occupied molecular orbitals of a pyridine radical. 


nated metal atoms, but the ligands break the octahe¬ 
dral symmetry. In complexes with D 3 symmetry, e.g., 
[Ru(bpy) 3 ]^+, the ligands have a reflection symmetry, cf. 
Fig. ^ ,d. For a single bpy ligand the highest energy 
ligand 7r-orbitals are even under this reflection whereas 
the lowest energy 7r*-orbitals are odd under the same 
reflection, as one would expect from simple symmetry 
arguments m- Therefore, linear combinations the tt- 
orbitals transform corresponding to the t 2 g representa¬ 
tion of Oh, whereas the 7r*-orbitals form a representation 
of tiu^ Therefore, 7 r-orbitals mix effectively with the oc¬ 
cupied metal d {t 2 g) orbitals but tt*- orbitals do not. Thus 
the highest occupied molecular orbitals (HOMOs), 
of the complex will have a significant contribution from 
both the ligand 7r-orbitals, and the metal-t 2 ^ or¬ 

bitals, d'^ . Neglecting smaller contributions from other 
ligand or metal orbitals, we have 

K, - cos 0 + TTt™ sin 6>, (1) 

where 0 parameterises the degree of mixing, and m G 
{1,2,3} labels the ligands and symmetry equivalent lin¬ 
ear combinations of d-orbitals. In contrast, the lowest 
unoccupied molecular orbitals (LUMOs) of the complex 
will be almost pure ligand tt* -orbitals. 

Low energy excited states can be well approximated 
by a single hole in the HOMO manifold and a single elec¬ 
tron in the LUMO manifold [17]. As both the HOMOs 
and LUMOs of the complex are three-fold degenerate one 


can label such states by two I = 1 pseudo-angular mo¬ 
menta, which we denote Lh and Ll respectively. We will 
only discuss this assignment for three real space HOMO 
spin orbitals, K^g - it is trivial to extend the following 
analysis to the LUMOs. By referring to these states as 
‘HOMOs’ and ‘LUMOs’ we are adopting the language of 
molecular orbital theory. However, we note that so-long 
as the are three local states related by rotations of 
27r/3 the discussion below goes through regardless of the 
degree of correlations in the states. It is therefore conve¬ 
nient to work in second quantised notation, so we define 
(r|aj^|0) = , where |0) is the ground state, \r) is the 

state with a hole at position r; spin labels are supressed. 

We introduce three ‘Bloch’ operators defined by 


b+ = sgn’^{-k)^Yl 

V rnn 


^i27rkml3 


( 2 ) 


where k G { —1,0,1}. Finally we identify the states cre¬ 
ated by the Bloch operators with the eigenstates of 
i.e., {0\bkLfjbl.\0) = k. The phase pre-factors [sgn^(—/c)] 
in the definition of the Bloch operators are required to al¬ 
low this assignment and maintain the required behaviour 
under time reversal symmetry. 

As the LUMOs are pure ligand orbitals the exchange 
interaction will be dominated by the exchange interaction 
between the ligand tt and 7 r*-orbitals, In contrast the 
SOC on the metal, is much stronger than the SOC 
on the ligands. Therefore states with one hole in the 
HOMO and one electron in the LUMO are described by 
the Hamiltonian 


Ho = JSh • Sl + \Lh • Sh, (3) 

where Sh is the (net) spin of the electrons in the HOMO, 
is the spin of the electron in the LUMO, J sin^ 0 
and A cos^ 0. Thus we expect positive J and A. 

If the excited state is sufficiently long lived for the 
geometry to relax it will be unstable to a Jahn-Teller 
distortion, which lifts the degeneracy. In terms of the 
pseudo-angular momenta this can be represented via the 
terms 

Hjt = 5Q +7Q [(Llf - (Llf] 

(4) 

where Q is the coordinate of the rhombic distortion per¬ 
pendicular to the C 3 -axis of the complex, 6 ( 7 ) is the 
coupling constant to the HOMOs (LUMOs) and k is the 
spring constant of the Jahn-Teller mode. 

It is helpful to briefly discuss the Jahn-Teller effect in 
the pseudo-angular momentum language, as this is not 
entirely intuitive. Consider the term 

= + (5) 

In terms of the Bloch operators = 2{b\bo + ^ 0 ^- 1 ) 
and = 2{b^_^bo + ^ 0 ^ 1 ). Hence, 
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Eiti [cm“^] 

Eiijii [cm 

Tj (l/fc^j) [ns] 

Tii (i/fcj/) [ns] 

Tin [/is] 

Ir(biqa)3 

14 

64 

107 (114) 

5.6 (5.7) 

0.36 (0.38) 

Ir(ppy )3 (in PMMA) 

12.2 

113 

154 (175) 

15 (17) 

0.33 (0.34) 

Ir(ppy)3 (in CH 2 CI 2 ) 

19 

151 

116 

6.4 

0.2 

Ir(dm-2-piq)2 (acac) 

9.5-10 

140-150 

80-124 

6.5-8.6 

0.33-0.44 

[Os(phen)]2(dppm)]^^ 

16 

106 

95 

13 

0.6 

[Os(phen)2(dpae)]^^ 

21 

92 

100 

10 

0.7 

Ir(piq)(ppy )2 

16 

91 

64 

10.5 

0.3 

Ir (4,6-dF ppy) 2 (acac) 

16 

93 

44 

9 

0.4 

Ir(pbt) 2 (acac) 

6 

97 

82 

25 

0.4 

Ir(piq) 2 (acac) 

9 

87 

47 

8 

0.3 

[Os(dpphen) ^ (dpae)] 

19 

75 

92 

9 

0.7 

[Os(phen)2(DPEphos)]^^ 

16 

68 

104 

14 

0.9 

[Os(phen)2(dppe)]^^ 

19 

55 

107 

12 

0.9 

R(piq) 2 (ppy) 

9 

56 

60 

6.4 

0.44 

[Os(phen)2(dppene)]^^ 

18 

46 

108 

15 

1.1 

[Ru(ppy)3]^+ 

8.7 

52 

230 

8 

0.9 

Ir(piq)3 

11 

53 

57 

5.3 

0.42 

Ir(4,6-dFppy)2(pic) 

9 

67 

47 

21 

0.3 

Ir(thpy) 2 (acac) 

3.5 

31 

113 

35 

1.5 

R(ppy)2(ppy-NPH2) 

6 

21 

188 

19 

1.8 

Ir(ppy-NPH2)3 

6 

20 

177 

15 

1.4 

Ir(ppy)(ppy-NPH 2)2 

6 

17 

163 

20 

2 

Ir(btp) 2 (acac) 

2.9 

22 

150 

58 

2 

Ir(btp) 2 (acac) 

2.9 

11.9 

62 

19 

3 

Ir(sl-thpy )2 (acac) 

3 

13 

128 

62 

3 

Ir(ppy) 2 (CO)(Cl) 

< 1 

< 1 

300 

85 

9 

[Rh(ppy)3]®+ 

- 

- 

4.5 X 10^ 

1.35 X 10^ 

650 


TABLE I: Key spectroscopic data for pseudo-octahedral d®-complexes. Ejji is the energy gap between the two lowest energy 
substates of Ti, Euju is the energy gap between the two highest energy substates of Ti and the total lifetime of substate 
m Tm = where and k^R and the radiative and non-radiative lifetimes of substate m. For Ir(ppy )3 and 

Ir(biqa )3 we also list 1/k^ (in bold) which, unsurprisingly given the high photoluminescent quantum yields in these complexes, 
shows the same trend as Tm- We are not aware of measurements of k^ in other relevant complexes. Note that in all complexes 
Eiji < Eli jII and ri > rn > rm, which suggests that < k^ < k^R^. To avoid selection bias we have included all and 
only those pseudo-octahedral d®-complexes included in Table 2 of the recent review by Yersin et ah [T]. The two rows for 
Ir(btp) 2 (acac) correspond to different sites. 






= b\h-i + 


3 L' 


2a\ai — 020.2 — a 3«3 + 2 ( 0^3 + 03 ^ 2 ) “ o\o 2 — o\ai — ajas — o\ai 


(6) 


It is therefore clear that this physics of Hjt is that of 
the T X t Jahn-Teller problem [or, once trigonal terms 
are included, below, the {A ^ E) x e pseudo-Jahn-Teller 
problem] and that this distortion corresponds to the so- 
called Eq distortion in the notation of, e.g., section 5 of 
Ref. [18]. The E^ distortion corresponds to terms pro¬ 
portional to (1/2)where 
curly brackets indicate anticommutation. 


In general the Jahn-Teller distortion could also induce 
a trigonal component of the distortion [which would cou¬ 
ple to (I/^)^, where u = H or L], however this does not 
produce any qualitatively new features and so, for sim¬ 
plicity, we neglect it below. Thence, the form of Hjt is 
constrained to the form given above by symmetry as: (1) 
terms that are proportional to odd powers of where 
P = y or z break time reversal symmetry and so may 
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not appear in the Hamiltonian for scalar Q and (2) for 
I = 1 any even power of is proportional to . 

However, the complexes in Table [I| are not octahedral, 
but trigonal. In terms of the pseudo-angular momenta, 
this introduces the additional terms 

Ht = A{L^H?+T{Ll)\ ( 7 ) 

where A (T) is the energy differences between the HOMO 
and HOMO -1 (LUMO and LUMO+ 1 ) in the trigonal 
ground state, geometry. Indeed, it immediately fol¬ 
lows from time reversal symmetry that the trigonal terms 
in the Hamiltonian are constrained to take this form. 
t 2 g ai e and tiu a 2 + e on lowering the sym¬ 
metry from Oh to D 3 . Therefore, the two pairs of e 
states are allowed to weakly mix, stabilising the = ±1 
states and destabilising the = ±1 states. Thus one 
expects that both A and T will be positive [13]. The 
approximate D 3 symmetry of the complexes with lower 
symmetry, e.g. C 3 , complexes considered here means 
that we expect both parameters to remain positive for 
all of the complexes considered here [19]. Thus the ef¬ 
fective pseudo-angular momentum Hamiltonian for the 
low-energy excitations is 

H = + (8) 

By definition Q = 0 in the Sq geometry and, by suit¬ 
ably rescaling the parameters, one may define Q = 1 in 
the Ti geometry. Similarly any trigonal component to 
the Jahn-Teller distortion can be taken simply to shift 
the value of A (T). Therefore, up to constants, in the Ti 
geometry the effective electronic Hamiltonian is 

H = JSH-SL+XLH-SH + A{LI,f+r{Llf 

+5[iL-Hf - + liiLlf - (Lin (9) 

As well as describing systems displaying a Jahn-Teller 
distortion, this model is also appropriate for heterolep- 
tic complexes. Indeed for appropriate choices of A, J, T 
and 7 one can parameterise arbitrary energy differences 
of the frontier orbitals. We discuss the values of these 
parameters in the Appendix. On the basis of this dis¬ 
cussion, for Ir(ppy) 3 , we take A/J = 0.2 and A/J = 0.5, 
with J ^ 1 eV; J < A and 7 ^ T below. Clearly, for 
example, A is strongly dependent on the transition metal 
in question. However, our main qualitative results are 
insensitive to the values of these parameters - to empha¬ 
size this we explore a wide range of other parameters in 
the sup. info. 


III. RESULTS 
A. Octahedral model 

Before considering the full pseudo-angular momentum 
model, iJ, it is important to understand the symmetries 



A/J 

FIG. 2: Energy eigenvalues of Ho for states with = 0. 
At A = 0 the singlets have E — 3J/4 and the triplets have 
E = —JjA. For A > 0 the labels “singlet” and “triplet” are no 
longer strictly defined (in their usual sense) nevertheless the 
relatively small energy shifts suggest that these labels retain 
some meaning, this claim is supported by directly examining 
the character of the eigenstates. It is interesting to note that, 
already in the octahedral problem, the lowest energy (non¬ 
degenerate) state has no singlet contribution to its wavefunc- 
tion for any value of A, thus radiative transitions from this 
state are forbidden. 


of Hq. (i) Ll does not couple to any of the other vari¬ 
ables. Therefore, L\ and are good quantum num¬ 
bers. (ii) We can define a ‘total’ angular momentum, 
I = Lh S', where S = Sh E Sl. P and com¬ 
mute with Ho therefore / and H are also good quantum 
numbers. 

We plot the energies of the exact solutions of Hq in 
Fig. (Table [H] gives the basis used for all calculations 
in the paper). For simplicity Fig. shows only the so¬ 
lutions with = 0 - because is decoupled from the 
other angular momenta it can be immediately seen that 
the other solutions simply triple the degeneracies of all 
states. Note that, firstly, the spectrum of Hq is not very 
similar to those of the pseudo-octahedral complexes we 
are seeking to model. However, this model is an impor¬ 
tant stepping stone to understanding the full Hamilto¬ 
nian. Secondly, the eigenstates can be classified by their 
total angular momentum quantum number, /, and, as 
Ho is SU(2) symmetric, have the expected 2/ -j- 1 degen¬ 
eracy. Thirdly, all of the singlets have / = 1; as Lh = 1 
and, by definition, singlets have S' = 0. This means that, 
regardless of how strong the SOC is, the singlets can only 
mix with the 1 = 1 triplets. Therefore radiative decay 
from the / = 0 and 1 = 2 triplets is forbidden by the 
conservation of /. 


B. Trigonal model 

In Fig. I^we plot the spectrum of the trigonal model 
with no Jahn-Teller distortion. Ho E Hf. Again, for sim¬ 
plicity, we only show the solutions with = 0. In this 
case each state has partners with = ±1 that have en- 
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Name 

Relationship to eigenstate 
of Ho when A = 0, 

r 

X^ 

\StStiLj,) 

Singlets 
mixed with 

in full model 


|1,0,0) 

-1 

1 

^ (1 -1 ;t^» 


- 

\ Sx ) 

^(|1,1,0)-|1,-1,0)) 

1 

-1 

1 (1 m) - 1 nt) -1 +1 w)) 


- 

\Sy) 

;T(|l,l,0) + |l,-l,0)) 

-1 

-1 

1 (1 t;t) -1 m) +1 tw) -1 w)) 


- 

\Tl) 

|0,0,1) 

1 

1 

(1 W> + 1 nt» - (1 + 1 It^)) 

None 

\T.) 

|1,0,1) 

-1 

1 

(1 tw -1 nt» 


15.) 

\n} 

^(ll,l,l)-ll,-l,l)) 

-1 

-1 

1 [i n^> -1 - 75 (1 nt) +1 ;n) -1 tw -1 

w»] 

\^y) 

\Ty) 

^(ll,l,l) + ll,-l,l)) 

1 

-1 

5 [i tt^> +1 (1 nt) +1 ;tt> +1 m) + \ 

w»] 

\S.) 

K-} 

|2,0,1) 

1 

1 

(1 w> +1 nt» + (1 +1 


None 

IT,,) 

^(|2,1,1)-|2,-1,1)) 

-1 

-1 

5 [i n^> -1 ^ (1 nt) +1 ;n> -1 tw -1 

w»] 

\^y) 

IT,.) 

^(|2,1,1) + |2,-1,1)) 

1 

-1 

5 [l tt=^) + 1 4-1'=^) + (1 tilr) + 1 itlT) + 1 ti-IJ-) + 1 

w»] 

\s.) 

|T,y) 

^(| 2 , 2 , 1 )-| 2 ,- 2 , 1 )) 

-1 

1 

(1 tn> -1 ;w) 


\s.) 

|T,2-,2) 

^(| 2 , 2 , 1 ) + | 2 ,- 2 , 1 )) 

1 

1 

(1 ttt> +1 


None 


TABLE II: The basis set used in this paper. The wavefunctions are given in the form \SIShLh) with t (i) indicating 
= +1/2 (—1/2) and ^ and indicating = 1,0 and —1 respectively. = (~ 1 )^ • fhis table we list only the 
= 0 [C^ = (— 1 )^^ = 1 ] states. Each state has two partners with = 1 and hence — — 1 . The latter two, but not the 
former, mix under the action of the full Hamiltonian. 


ergies that are higher by T and display twice the degen¬ 
eracy of the = 0 state. The trigonal terms break the 
SU(2) symmetry of the octahedral model and therefore 
lift the three- and five-fold degeneracies. The calculated 
spectra are now like those calculated from first-principles 
for relevant complexes. For example, if trigonal symme¬ 
try is enforced for, e.g., [Os(bpy) 3 ]^+, Ir(ppy) 3 , Ir(ptz )3 
relativistic TDDFT calculations predict that SOC splits 
Ti into a non-degenerate state (I) and, at slightly higher 
energies, a pair of degenerate states (II and III) [811^ITT]. 

We saw above that in the octahedral model radiative 
decay from the lowest energy excited state (I-^O) is for¬ 
bidden by the conservation of I. Because Hf breaks the 
SU(2) symmetry of the octahedral model P no longer 
commutes with iT, nevertheless and remains a 
good quantum numbers for the trigonal model. Fur¬ 
thermore, the Hamiltonian is time reversal symmetric, 
therefore the parity of an eigenstate under time reversal, 
T = ±1, is also a good quantum number. Note however, 
that P does not commute with time reversal so it is not, 
in general, possible to form states that are simultaneously 
eigenstates of both. However, one may define states that 
are simultaneous eigenstates of the T and = (—1)^ . 
Therefore, we take these as our quantum numbers, cf. 

Table m 

For all parameters studied substate I is composed of 
the basis state |Ti) admixed with |T^ 2 ) and has quantum 
numbers = T = +l, I/|^ = 0 whereas states II and 
III are a degenerate pair with = —1, T = ±1, = 

0, cf. Fig. whose largest contributions come from 
\Tx) and \Ty). The singlet states with the same quantum 


numbers contribute to substates II and III, but all of 
the singlets are forbidden from mixing with substate I 
by the combination of time reversal symmetry and the 
conservation of Hence, the 1^0 transition remains 
forbidden in the trigonal model. Both experiments pi [20] 
and relativistic TDDFT calculations isti] find that the 
radiative rates for the transitions II—and III—are 
more than an order of magnitude faster than that for 
1^0, cf. Table [n] The small non-zero decay rate for 1^0 
may arise from either Herzberg-Teller coupling mm or 
mixing of state I with higher energy singlet states, which 
are not included in the pseudo-angular momentum model 

EHII]. 


C. Full model 

Finally, we turn to the full pseudo-angular momentum 
model, H [Eq. ^]. does not commute with Hjt- 
However, ~ where the 

ladder operators are given by = Lfj ± there¬ 
fore P is conserved modulo two. Thus, X^ is conserved 
even for a trigonal system that has undergone a Jahn- 
Teller distortion, cf. Table [n] Similarly is conserved 
modulo two, which gives rise to the quantum number 

We plot the spectrum of = 1 states in Fig. 1^. In 
Fig|^, we plot the same results, but only show the three 
lowest energy substates, I-III, which are of primary tech¬ 
nological interest. One sees that the although states II 
and III are degenerate at = 0, a Jahn-Teller distortion 
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FIG. 3: Solution of the pseudo-angular momentum model 
of a trigonal complex, a) spectra for A = J/5 and varying 
A/J; b) spectra for A = J/2 and varying A/J. The quantum 
numbers of the states are also indicated. In both panels the 
states with quantum numbers labelled as T = ±1 are two-fold 
degenerate. The eigenstates with = ±1 (not shown for 
clarity) have the same properties except that their energies 
are increased by T and all of the degeneracies are doubled 
corresponding to the two values of = ±1. 


rapidly lifts this degeneracy and for reasonable values of S 
one finds that there is a much smaller energy gap between 
substates I and II than between II and III. This is what 
is observed experimentally [DEolIll] in a huge range of 
complexes (Table IS. We will see below that this splitting 
is the signature of the localization of the excitation to a 
single ligand. 

Note, in particular, that substate I remains an admix¬ 
ture of |Ti) with |T^ 2 ) and and has quantum 

numbers T = = +1. As none of the singlet 

states have these quantum numbers, cf. Table |n| this 
state is forbidden from mixing with any of the singlet 
states in the model by conservation of T, and . 
Therefore substate I remains a pure triplet and is forbid¬ 
den from decaying radiatively, irrespective of the strength 
of the SOC. 

The radiative rate of the mth eigenstates of the full 


FIG. 4: Solution of the pseudo-angular momentum model 
of a complex with broken trigonal symmetry - due either to 
chemical modihcation or excited state localization. Panel (a) 
shows the full spectrum for states with Panel (b) 

shows only the Ti substates, which are our primary concern. 
Here we take A = J/2 and A = J/5. 


Hamiltonian, is given by 


E {SMSn){Sn\^m) (10) 

/3 G {x,y,z} 
n e {x,y,z} 


where is the excitation energy of the mth state, and 
a is the fine structure constant. Because of the underly¬ 
ing octahedral symmetry we take {So\fi^\Sn) to be inde¬ 
pendent of n and further we assume that the zero field 
splitting is small compared to the Sq Ti excitation 
energy, i.e., that Ej Ejj Em. It is also convenient 
to define 


g ^ Aniee'^aE^ 

^ 9(47reo)2ft3 


E (5o|Ai/5|5„) 

/3 G {x,y,z} 
n G {a:, y, z} 


( 11 ) 


this corresponds to the radiative decay rate for a pure 
singlet with an excitation energy equal to that of the Ti 
manifold. 

We plot the radiative decay rate in Fig. State I is 
dark - as expected from the conservation laws derived 
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FIG. 5: The radiative decay rates of the three substates of 
Ti. The conservation of T, and Cf' leads to the absence 
of radiative decay for state 1. It can be seen that once the 
Jahn-Teller distortion becomes significant the radiative decay 
rate from state II is significantly smaller than that from state 
III, in good agreement with experiment (cf. Table |^. Here, 
as above, we take A = J/2 and A = J/5. 



0.0 0.2 0.4 0.6 0.8 1.0 


^/A 

FIG. 6: The degree of localization in the three 

substates of Ti, v G {I, II, III}, where = 

I^L^Oct - I (a\^aia + al^a 2 a^ | Here, as above, 
we take A = Jl2 and A = J/5. 


above. Furthermore, once the Jahn-Teller distortion be¬ 
comes significant one finds that the radiative decay from 
state II is significantly slower than the radiative decay 
from state III. This is in precisely what is observed in ex¬ 
periments jpunun] on pseudo-octahedral complexes 
(cf. Table 0. 

It is straightforward to understand both the changes 
in energy and the radiative rates of states II and III. For 
J > 0 (J < 0 simple reverses these effects) the trigonal 
perturbation lowers the energy of (stabilises) states that 
are antibonding between the X^ = ±1 orbitals, e.g., \Sx) 
and \Tx), and raises the energy of (destabilises) those that 
are bonding between the X^ = ±1 orbitals, e.g., \Sy) and 
\Ty). It is clear from Table [ll| that whereas \Sx) and \Ty) 
are even under time reversalJXc) and |5'y) are odd. Thus, 
SOC mixes \Sx) with \Ty) and \Sy) with \Tx). Hence 
the trigonal distortion increases the energy difference be¬ 
tween the triplet and singlet basis states that contribute 
to state II (i.e., \Sx) and \Ty) for 6 > 0); whereas trigo¬ 
nal symmetry reduces the energy difference between the 
triplet and singlet basis states that contribute to state 
III {\Sy) and \Tx) for 6 > 0). Thus the symmetry of 
the model dictates that < kj^^, as is observed 

experimentally [D [201 m] , see Table |I| 

Finally, we turn to the question of localization in the 
excited state. To measure this we define 





- + al^a2a 


)|«). 


( 12 ) 


We plot for the three substates of Ti in Fig|^ The 
lowest energy excitation, I, is completley delocalized for 
J = 0 but rapidly localizes for J > 0. It is interesting 
to note that both Sp and Spi are non-zero for J = 0. 
However, for J = 0 states II and III are degenerate and 
Ell = consistent with trigonal symmetry. (A neg¬ 

ative value of E^ indicates that the state disfavors occu¬ 


pation of site 0.) Nevertheless, for J > 0, one observes 
a rapid increase in Eu whereas Em grows only rather 
slowly. 

It is therefore clear that the pseudo-angular momen¬ 
tum model predicts significant localization for values of S 
compatible with the observed experimental results that 
k]i < k^^ < k]^^ and Ei^u < X^ipm, cf. Table We there¬ 
fore conclude that all of the complexes in Table [T] show 
significant localization in their excited states. 


IV. CONCLUSIONS 

The pseudo-angular momentum model gives a natu¬ 
ral explanation of the zero-field splitting observed in a 
wide range of pseudo-octahedral organometallic com¬ 
plexes. Furthermore, the conservation laws, and hence 
selection rules, inherent in the model give a natural ex¬ 
planation of the relative radiative decay rates of the three 
sublevels of Ti. We stress that none of the results derived 
here rely on perturbation theory - therefore these con¬ 
clusions hold even when the departures from octahedral 
or trigonal symmetry are large. This immediately ex¬ 
plains why the properties of the Ti states are so similar 
in both homoleptic and heteroleptic complexes. Further¬ 
more, for parameters compatible with the experimentally 
measured energies and radiative rates of the substates of 
Ti, the pseudo-angular momentum model predicts that 
exciations I and II are strongly localised - although III 
remians well delocalised. Thus we conclude that all of 
the complexes in Table [T] show significant localization in 
their two lowest energy excited (sub)states. 

It is interesting to note that when the radiative rates 
of individual sublevels, have been measured, rather 
than excited state lifetimes, r^, it is found that the rel¬ 
ative rates are in good accord [HEO], cf. Table 1^ This is 
consistent with the high photoluminscent quantum yields 
observed in these complexes. This suggests that the non- 








radiative decay rates of the individual sublevels are deter¬ 
mined by similar conservation laws. Therefore, it would 
be interesting to investigate non-radiative decay rates in 
a suitable extension of the pseudo-angular momentum 
model. 

We note that the pseudo-angular momentum model de¬ 
scribed above can be naturally extended to understand 
the properties other molecules and complexes where the 
low-energy excited states correspond to transition be¬ 
tween degenerate or approximately degenerate states. 
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V. APPENDIX: ESTIMATION OF 
PARAMETERS 

While we will not make a detailed parameterization of 
this model - an idea of the relevant parameter ranges 
can be obtained from previous experiments and density 
functional calculations. Ir(ppy )3 has been particularly 
widely studied and so is an ideal material to compare 


with. 0 7 r /4 as the HOMOs are found to have about 

50 % metallic weight [H [TOl [TTl [ 22 ]. Nozaki [ 8 ] found 
that for Ir = 550 meV and sin^ 0 = 0.4, yielding 
A = 220 meV. Smith et al. [TT] considered the C 3 Sq 
geometry and found from ground state calculations show 
that the gap between the HOMO and HOMO -1 is 140 
meV (LUMO and LUMO+1 is 90 meV), which may be 
taken as an estimate of A (F). However several authors 
Nozaki [S1II71I13] have noted that the values of A and 
F are difficult to calculate from first principles - there¬ 
fore these numbers should be treated with some caution 
and are likely to be underestimates as interactions sig¬ 
nificantly increase the effective values of A and F. For 
ppy it has been estimated [23] based on the absorption 
spectra, emission spectra, and emission lifetimes [24] that 
Jtt ^ 2 eV; and for an isolated Ir ion ^ 0.43 eV [25] . 
Taking 0 ^ it/A yields J ^ 1.4 eV and A ^ 300 meV. For 
concreteness we take A = J/5 and A = J /2 in the main 
text. However, our results are insensitive to the values 
of these parameters - to demonstrate this we explore a 
range of other parameters in the sup. info. 6 and 7 are 
not straightforward to estimate from previous work and 
will be left as free parameters, however as the C 3 sym¬ 
metry remains evident even in the Ti geometry of the 
excited this suggest that 6 (7) is not significantly larger 
than A (F). 
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